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ABSTRACT 
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2-jet cross sections are considered as a background to 77 —>■ H* —>■ bb (standard model). 
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and t-quark production, the higher order QCD corrections are found to be significant. 
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I. INTRODUCTION 



Higher (next-to-leading) order corrections (HOC) for heavy-quark {Q,Q) production 
in unpolarized particle collisions have been determined in detail. Bi For polarized particle 
collisions, however, analytical results were still absent. Even for the unpolarized case, only 
virtual -|- soft corrections have been presented analytically. i Apart from general reasons, 
well known from unpolarized reactions, knowledge of HOC for Q, Q production in polarized 
processes is important for several special reasons. 

Beginning with polarized 77 collisions, which is the subject of the present work, one 
reason of special interest is the following. A 77 collider becomes particularly important for 
searches of the standard model Higgs boson when its mass is below the W~^W~ threshold. 
Then the predominant decay is H ^ bb and the background comes from 77 — > bb with direct 
or resolved photons. Leaving aside the latter, for the moment, use of polarized photons of 
equal helicity (when the angular momentum has Jz = 0) suppresses this background by 
a factor m^/s. Si This holds, however, only for the lowest order of as- HOC necessarily 
involve the subprocess 77 bbg, and gluon emission permits the bb system to have J 7^ 
without suppression; this may result in a sizable background. Of course, another reason 
the Jz = channel is important is that the Higgs signal comes entirely from it. Thus, we 
maximize the Higgs to background ratio in two different ways. 

Furthermore, at higher energies, it will be possible to produce top-quarks in photon- 
photon collisions. This, when combined with other data on top-quark production from e~^e~ 
and pp collisions, should certainly improve our knowledge of the top-quark parameters. The 
HOC could have a significant effect on the threshold behaviour. It is also interesting to 
examine the spin dependence of the HOC in this region. 

In this paper we present complete analytical results for heavy quark production by both 
polarized and unpolarized photons. Numerical results are presented for 2-, 3- and 2+3-jet 
cross sections for the cases where the initial photons have total spin Jz = and Jz = ±2. 
For 6-quark production, this is analyzed as a background to Higgs production. We also 
consider t-quark production for energies not too far above threshold. 

The analytical results presented here are also useful in determining the production 
of heavy quarks in polarized photon-proton (proton-proton) collisions. This is because 
the process 77 QQ{g) is the Abehan (QED) part of the subprocess •yg QQ{g) 
{99 ~^ QQi9))i which is by far the dominant subprocess in 7-p (p-p) collisions. i'i The 
non-Abelian part of ^g — > QQ{g) (gg — > QQ{g)) remains to be calculated. 

II. LEADING ORDER CROSS SECTIONS 

The contributing graphs are shown in Fig. 1. We introduce the variables (momenta as 
in figure) 

s = (pi+p2)^, t = T - = {pi - p^y - m^, u = U-m^ = {p2-p:iY-w? (1) 
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Figure 1: Lowest order contributions to 77 QQ. 

and 

S2 = S2 - = {pi +P2- Pzf -m^ = s + t + u. (2) 

where m is the heavy-quark mass. Defining 

t —u , , 

+ u,.— (3) 

we may express 

t = —s{\ — v\ u = —svw, S2 = sv{l — w). (4) 
The polarized and unpolarized squared amphtudes are defined, respectively, as 

A|Mp = i(|M(+, +)p - |M(+, |Mp ^ 1(|M(+, +)p + |M(+, (5) 

where M(Ai,A2) denotes the Feynman amplitude with photons pi, p2 having helicity Ai, 
A2 respectively. The same holds for the cross sections. 

For the examples considered in this paper, it is of interest to calculate (numerically) the 
cross sections for a specific helicity state, (t(Ai,A2). We present analytical results for the 
polarized and unpolarized cross sections Act, a. From we can obtain the desired cross 
sections via 

(t(+, +) = (T + Ad, -) = a - Aa. (6) 

Defining 

we may express the n (= 4 — 2£:)-dimensional 2-body phase space as 

[A]^ = K{e) [(2m)^[A]|M|L2] '^(1 - (8) 

It will become necessary to work in n dimensions when we determine the HOC (see next 
section for details). 

The resulting leading-order (LO) cross sections are, in n dimensions, 

avdw ^ \ tu tu \ tu \ 
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Figure 2: Loop graphs for 77 — > QQ. (a)-(c) self-energy diagrams; (a')-(c') mass counterterm 
diagrams corresponding to the graphs (a)-(c); (d),(e) vertex diagrams; (d'),(e') dimensional 
reduction counterterm diagrams corresponding to graphs (d),(e); (f) box diagram. 



where Nc (=3) is the number of quark colors and eg is the fractional charge of the heavy 
quark. Making use of (H), (||) we see explicity that (i(TLo(+, +)/dvdw is suppressed by order 



III. LOOP CONTRIBUTIONS 



The loop contributions arise from the diagrams of Fig. 2 and their pi ^ p2 interchange. 
These diagrams contain both ultraviolet and infrared singularities. To regularize them, we 
use dimensional reduction!, where the momenta are in n dimensions and everything else 
is in four dimensions. This facilitates the handling of the Levi-Civita tensor et^'^^p. As we 
will show below, the analytical expressions for the cross sections are regularization scheme 
independent once all the contributions (including the gluonic bremsstrahlung) are added. 
Throughout, we work in the Feynman gauge. 

The heavy-quark mass and wave function renormalizations are performed on-shell. The 
self-energy graphs are shown in Figs. 2 (a)-(c) and the corresponding mass counterterm 
diagrams in Figs. 2 (a')-(c'). The factor 1/2 multiplying (b)-(c') comes from wave func- 
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tion renormalization. The bare mass and wave function are determined in terms of the 
renormahzed ones via 

mo = Zmtrir, ^0 = ^2^^^r, (10) 

where and Z2 are the mass and wave function renormahzation constants. Define 

where fi is an arbitrary mass scale which enters via the couphng in n dimensions: g —>■ gix^. 
In dimensional reduction we find, to order (7^, 

Z^ = l- ?>g^CeCF + ^) , Z2 = l- g'^CeCp + 5 + ^) (12) 

with Cp = 4/3. We use 1/e' to indicate which terms are of ultraviolet origin. 

In dimensional reduction we must add to the vertex diagrams of Figs. 2 (d), (e) appro- 
priate counterterms (d'), (e') in order to satisfy the Ward identity^ 

Zi = Z2 (13) 

between vertex and self-energy graphs, with Zi denoting the vertex renormalization con- 
stant. The Feynman rule for this vertex counterterm is found to be (in n dimensions) 

-q^ 1 

7'' — —Cp-f^- (14) 

with 

7^ = ig'" - Ca"-,. = C9,u = n; (15) 

here g^'^ represents the n-dimensional metric tensor with, formally, n < 4. 

When all the contributions to the physical cross section (including gluonic bremsstrah- 
lung) are added, the result is free of infrared divergences as there are no collinear singulari- 
ties here. Thus the only scheme dependent part might come from the vertex and self-energy 
graphs. Having satisfied the Ward identity (|13[) though, means that the scheme dependent 
part of the corrections cancels between vertex and self-energy graphs. This was explicitly 
verified by calculating the vertex and self-energy graphs in dimensional regularization. We 
also checked explicitly that there are no differences between reduction and regularization 
arising from any other contributions. More specifically, to obtain the dimensional regular- 
ization result for any particular contribution given in this paper, simply replace the LO 
term by the corresponding LO term from dimensional regularization. When all the contri- 
butions are added, the scheme dependent part of the LO term cancels along with the 1/e 
infrared divergence multiplying it. Hence, the absence of collinear divergences or vacuum 
polarization graphs leads to scheme independence. 

As was stated in Refs. ^ and a counterterm like (|l^) was used to remove an unphys- 
ical term. General one-loop counterterms have been developed^ to convert unpolarized 
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dimensional reduction results into the corresponding dimensional regularization results for 
the purely massless case. Also, certain equivalences between dimensional reduction and 
dimensional regularization have been noted. In the present case however, satisfaction of 
( [I3D is sufficient to ensure scheme independence. 

Adding the contributions of Figs. 2 (a)-(e') (and the t ^ u interchange) resulted in the 
ultraviolet finite vertex plus self-energy cross section 

= vra.CFC,— — + (5(1 - w) {2Ai 4(^(2) - Li2(— ))(1 + 3—) (16) 

avaw e avaw (47rj^ t 

- ln(- A)(8 A^H-—^) + A,{U,{—^) - e(2)) + A, + {t^ u)} 
where 

C = {ATryCpNcasa'^e^Qfi^'. (17) 

The corresponding polarized cross section, Ada^sc/dvdw can be obtained by replacing the 
Ai and dai,o / dvdw in (p!6| ) by AAi and Adaio/dvdw, respectively. The [A]Ai are given 
in Appendix B. We will use this notation throughout. We note the term ~ 1/e in ([16D 
representing an infrared divergence. Also, note that [A]Ai is proportional to the LO squared 
amplitude without the t ^ u interchange (see Appendix B, Eq. p^) . 

Since [A]dai^o/dvdw is in general regularization scheme dependent to 0{e) (working in 
n dimensions), we see explicitly that truly scheme independent cross sections will result 
only when all contributions are added and all infrared divergences are cancelled. 

In order to evaluate the box graph of Fig. 2 (f), we must reduce the resulting tensor 
integrals to scalar ones (conveniently listed in Ref. ^|) using projective tensor techniques. 
The tensor integrals have the general form 

D^'^^'^^'^'^'^\k,,k2,h,rm,m2,m3,m,) = (18) 
d'^q l,g^,gV,gVg^ 



{2%)^ {q - mi)2[(g + k^y - ml][{q + ki + kiY - ml][{q + k^ + k2 + h 



where the ki are general momenta. As an example, the vector box integral we encounter 
has the decomposition 

D^iPi, -P2, -Pi, 0, m, m, m) = p^Du - P2D12 - Pi -Dis- (19) 

In general, the scalar coefficients Dij are not independent. This simplifies somewhat the 
calculation. Noting that 

D^{p4, -p2, -pi, 0, m, m, m) = -D''{p3, -pi, -p2, 0, m, m, m), (20) 

we obtain 

D,2 = 1^11 - /^13, (21) 

since the Dij in both integrals are the same, due to the fact that they are scalars. Using 
the same approach, we reduce the number of independent D^j from seven to five in D'^'^ 
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and from thirteen to eight in D^"^. This method was quite helpful in keeping the very 
large intermediate expressions as short as possible. 

Adding the contribution of Fig. 2 (f) and the t ^ u interchange gives the virtual box 
cross section 

^ = 167ra,C^C',^^^^{21n(x)[l-ln(/5)]+2Li2(-x)-2Li2(a;)-3e(2)} 
dvdw dvdw sp 2e 



+ ^(1 - w) ^/^',7 {-8g r"" ^ " ln(x)ln(-t/m^) +2^[ln(x)(41n(l + x) -ln(x) 
[Arr)'' sp p 

- 41n(-t/m2)) +4Li2(-x) +2^(2)] + 253ln^(x) +4-^ln(x) +455ln(-t/m2) 

P 

+ 8BGU2{T/m^)+4Bj^{2)+4Bs + {t^u)} (22) 



where 

e(2) = y, P^^l- 4m7., (23) 

The [A]Bi are given in Appendix B. We see again the infrared divergence ~ 1/e. 

Independent calculations were performed using FORM@ and REDUCE0. The latter 
proved useful in factoring the expressions and cancelling powers in the denominators. 



IV. GLUONIC BREMSSTRAHLUNG CONTRIBUTIONS 

The bremsstrahlung diagrams are shown in Fig. 3. Squaring these diagrams (plus their 
Pi ^ P2 interchange), we obtain the 2 — 3 particle squared amplitude 



{2mf\M\ 



2-*3 



ei , ; , ^3 , 64 , P2-k ^ eg ,67 



C[-^ + 62 P2 ■ k ^ \ 2 + ^5 2 r + 



si P2 -Pi P2- Pi Pi - pi S2 P3- k P3 - k^ 

68 , P2-k'^ P2-k en/s2 , eu 

+ +69 + eio + - + 

Pi ■P4P2- Pi Pi - Pa Pi - Pa P2 ■ Pa Ps ■ k P2 ■ Pi Ps ■ k ^ 

, 614 , P2-k^ P2-k 

+ 2 r ^ 2 + ^15 r + — r + (pi ^P2,t^ u) (24) 

P2 -Pi P3-k P2-Pi P3-k^ P3-k ps-k 

As before, we may obtain A|M|2^3 by replacing the Cj in (^) by Ae,. The [A]ej are given in 
Appendix B. Again, independent calculations were performed using FORM and REDUCE. 
The former proved useful in partial fractioning and other reductions of the dot-products. 

To obtain the total bremsstrahlung contribution to [A]da/dvdw, we perform the phase- 
space integrals in the frame where pi and k are back-to-back. We find (in agreement with 
Ref. H) for the 2—^3 phase space 
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where 



and 



^ 1 



w) 



l-2e 



J Jo Jo 



(26) 
(27) 



The gluon angles 6i and 62 are defined in Appendix A along with all the momenta. 

We first evaluate all the phase space integrals in four dimensions since, for w 7^ 1, all the 
integrals are finite. For w = 1, the terms in (^4]) with coefficients Cj are singular through 
the relationlll 

(1 - w)-'-'^ = -1^(1 -w) + — ^ + Oie) (28) 
2e (1 — w)+ 

where the function 1/(1 — 1;;)+ is defined through 

/H-/(i) 



1 f{w) 

dw- r 

wi [1 — W) 



dw- 



Wl 



W] 



+ /(l)ln(l-u;i) 



(29) 



This means that, for these terms, the integrals must also be evaluated in n dimensions in 
the limit w ^ 1, keeping their 0{e) part. The resulting integrals are straightforward. 
The final result is, with y = \J {t + uY — Am'^s, 



da 



Br 



dvdw 



2S2 



■In — 63 + 



4^2 



(47r)2 45*2 S2{s + u) w? rn?{s + uy 

+ 65/5 + es/g + 69/9 + eio/10 + 613/13 + 615/15 + ei6/i6 + {t^ u)} 
1 

\p, A In Pn -J- ZL — : 

da-Lo 1 



64 



+ 



CK(0) 1 2S2, T + U-y S2 

n ^ — TT^^i^i + — In Cq + 4— 67 + ^2/11611 

[l — w)^ (47r)^ S2 y T + U + y 



+ ^2/12612 + 52/14614 + (t ^ n)} + ^TxasCpC, 



X 



21na;(21n 



sv 



e' 



1 - 2 Li2 



where the integrals /j are given in Appendix C. 



dvdw sf3 
4/9 \ ^, 2 ' 



{(2m^ 



+ 2s(3 



l-21n^ + l 



(30) 



V. PHYSICAL CROSS SECTIONS 



We may obtain the 2+3-jet cross section by adding (H), (|T6D, (^2]) , and ( ^OD 

,da2+z 



[A] 



dvdw 



dvdw 



dvdw 
7 



+ [A]— ^ + [A]- 



dvdw 



dvdw 



(31) 




Figure 3: Gluonic Bremsstrahlung graphs for 77 QQg- 

We notice the cancellation of all the 1/e infrared divergences, leading to a finite, scheme 
independent result. 

At this point it is useful to note that for s S> 4m^, the LO cross sections are large in 
the forward and backward directions. Since jets going down the beam pipe are difficult to 
measure experimentally, angular cuts are necessary for bb production well above threshold. 
At the same time, we reduce the bb background to the Higgs signal. This also helps eliminate 
resolved photon contributions where the partons within the photon participate, as opposed 
to the direct contributions, which we present, where the photon is structureless. This is 
discussed at the end of this section. 

Let 63 denote the angle between p^, and pi in the 77 cm. Then the integrated 2+3-jet 
cross section, with the constraint | cos 6*3! < cos^c) ior some 9c, is given by 

[A](T2+3(s) = Hdv Cdw e(cos'0, - cos2^3)[A]^^ (32) 
Jvi Jwi dvdw 

where 

1 1 1 

2 2 s v{l — V) 

and 

COS 6*3 = ^ ' =. (34) 

w(l — f + vwY — Am? / s 

Alternatively, we may convert to da/dcosOsdw and integrate directly over 6^ and w. 
The integrated 3-jet cross section is given by 



[A]a3(s) = — \| Hdv f dw e{cosHc-cosH3)f{0) fdn(2m?[A]\M\ 
(47r)^ Jvi Jwi J 

X e((p3 + kf - ?/cutS)e((p4 + k)^ - VcutS) 



2^3 



(47r) 



'2 /"i(^2 ~ r 

2 I dv dw Q{cos^9c-cos'^e3)f{0) dn{2mf[A]\M\l^2 

vi Jw\ J 

\2 „^. ... _ 1 (2/cut - rr?/s) 



X e((p3 + fc) -1/cuts); W2 = l-^ (35) 

V 



The angular integral is given by (|27D with e = 0. The dot-products involved may be 
explicitly expressed as functions of f , w, 6i and 62 using the parametrizations of Appendix 
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A and Eqs. (|^). We have imposed the constraints, (pa + kY > ycnts and (p4 + k)"^ > ?/cutS. 
With a suitable choice of ?/cut we may simuhaneously cut out events with 2-jet topology 
and avoid the soft divergence. We effectively eliminate the soft and coUinear gluons from 
the 3-jet cross section, with the degree of softness and coUinearity being specified by ^cut-s. 
The desired 2-jet cross section is obtained by the difference 

[A]a2{s) = [A]a2+s{s) - [A]a,{s). (36) 

Since (J2+3 and (J3 are both infrared finite and separately observable quantities, this serves 
as a reliable and unambiguous method for defining a2- 

In discussing the numerical results, it will be convenient to split [A] 0-2+3 as follows, 

[A]a2+3 = [A]aLO + [A](Tg + [A](Th, (37) 

where [A] as represents the contribution to the HOC coming from terms proportional to 
(5(1 — w) and 1/(1 — w)+, and [AJcth represents the rest. In usual terminology, [A]crs 
represents virtual and soft contributions whereas [A](Th represents hard radiation. 

So far we have only considered direct contributions, i.e. no resolved photon contribu- 
tions. The reason is the following. Well above the QQ threshold, (T2+3 and 0-3 will certainly 
receive sizable resolved photon contributions. Now, resolved photon events are generally 
accompanied by a jet making small angles with respect to the beam axis. For the 2-jet 
cross section (which is of physical interest), experiment can reject resolved photon events 
(and other unwanted events) as being those for which the observed jets have total energy 
measurably lower than ^/s.u This is because, due to the angular cuts, experiment will not 
observe the jet making small angles. Hence, there will be missing energy. Of course, we are 
assuming a rather well defined initial photon energy, which may be experimentally difficult. 

For top-quark production, not too far above threshold, the resolved contributions will 
be negligible in all the cross sections. This is because the dominant resolved contribution 
comes from (77 — > QQ, where the gluon originates from one of the initial photons, having 
a fraction x of its momentum. Near threshold, the gluon will have to carry a large fraction 
of the photon's momentum; and for x ^ 1, the gluon distribution in the photon is highly 
suppressed. As well, 3-jet states arising from hard gluonic radiation will be suppressed due 
to the restricted phase space. The (near) absence of resolved contributions and the non- 
suppression of the Jz = cross section for 2 ^ 2 kinematics, not too far above threshold, 
implies that we needn't worry about whether the events are 2- or 3-jet (even though 3-jet 
events are either very seldom or none, depending on s). 

VI. NUMERICAL RESULTS 

Here we present numerical results for b- and t-quark production in next-to-leading order. 
Throughout, we evaluate asiQ"^) (2-loop) with Q"^ = s, A = 0.2 GeV and the number of 
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Figure 4: Cross sections for 77 bb{g): ulo (dotted line), 0-2+3 (dashed), (dash-dotted) 
and 0-2 (solid), with 6^ = 30° and ycut = 0.15 for 20 < ^ < 200 GeV; (a) a{+, +); (b) ai+, -). 

flavors taken as Np = 5 since we are well above the bb threshold. We take = 4.7 GeV 
and mt = 174 GeV.0 For 3-jet cross sections, we use ?/cut = 0.15. Some justification for 
this choice of ?/cut is in order. Experimentally, it is useful to have a small value of ?/cut so 
that for the 2-jet cross section we ehminate, as much as possible, events with 3-jet topology 
via (pBI). Theoretically, there are limitations. If one chooses ?/cut too small, then the 
infrared divergence ruins the perturbation expansion, since the 3-jet cross section becomes 
unphysically large. To control this, an all-orders resummation would be required. We find 
that ?/cut = 0.15 is the most suitable choice in light of the above considerations. 

Fig. 4(a) presents crLo(+, +), o"2+3(+, +), o"3(+, +) and cr2(+, +) for 6-quark production 
in the range 20 < i/i < 200 GeV with 9c = 30°. As expected, the LO cross section is highly 
suppressed for large ^/s, but not the 3-jet. In fact c73(-|-,-|-) makes a sizable contribution 
to cr2+3(+,+). Hence o"2(+,+) gets somewhat suppressed relative to cr2+3(+,+). For 
20 ^ ^/s ^ 40 GeV the corrections (72+3 — (Tlo are seen to be slightly negative. 

Fig. 4(b) presents the same cross sections for = ±2, i.e. crLo(+5— ), <^2+3{+, —), 
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Figure 5: Same as Fig. 4, except with 9c = 45^- 



(T3(+, — ) and o"2(+, — ). The major difference is that o"lo(+, — ) and (T2+3(+, — ) suffer no 
suppression at large a/s. Hence the 3-jet contribution to cr2+3(+, — ) is not so significant 
and o"2(+, — ) remains large. We also notice, CT2+3 ^<Jlo throughout. 

Figs. 5 (a),(b) present the same quantities as in Figs. 4 (a),(b) except with 9c = 45°. 
The major difference is that the cross sections are smaller everywhere and cr2{+,+) is 
particularly suppressed for 30 < ^/s ^ 60 GeV. This reflects the fact that the 2-jet events 
tend to occur at smaller angles. 

An interesting feature of the HOC arises for both (J2+3 and Acr2+3. In both cases, as 
and (Th are much larger that ctlo, for s ^ Aw?. But they have opposite sign and are of 
almost equal magnitude, leading to large cancellations. In other words, the "virtual + soft" 
part conspires with the "hard" part to yield HOC which are under control. 

Fig. 6 presents the 2-jet background to the Higgs decay 77 H* — * bb. We have used 
the standard model Higgs cross section of Ref. ^ which takes 6c = 30° and an average value of 
< A1A2 >= 0.8. The photons are produced by laser backscattering off electrons (positrons) 
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Figure 6: Two-jet bb background to standard model Higgs decay: 77 — > H* bb (solid line), 
(Tlo (dotted) and (T2 (dashed) for 20 < ttt-h < 200 GeV. Number of Higgs events taken from Ref. 
1^. Here 9c = 30^, < A1A2 >= 0.8 . The other experimental parameters are described in the text. 
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Figure 7: Cross sections for 77 tt{g): ctlo(+) — ) (lower dashed line), (T2+3(+, — ) (lower solid), 
clo(+i+) (upper dashed) and o"2+3(+,+) (upper solid) for 1 < -y/i/2m < 1.4: (a) 6c = 0; (b) 
Be = 30°. 



at an e~^e~ collider with E^+e- = 500 GeV. As well, Ref. ^ uses an effective integrated 
luminosity of Lcs = 20 fb~^ and a 77 energy spread of Fexpt = 5 GeV; y/s = mn ± rcxpt/2. 
Using the expression of Ref. ^ for converting the 77 — > bb{g) cross section into number of 
events, we obtain the LO and 2-jet curves shown in Fig. 6. 

At large y/s, the increase in (T2{+, +) relative to o"lo(+5 +) is compensated by a decrease 
in <72{+, — ) relative to crLo(+, — ), so that a2{< A1A2 >= 0.8) doesn't change radically. In 
the end, the 2-jet cross section is still well below the Higgs signal for 90 ^ mn ^ 150 GeV. 
With higher degrees of polarization, we could do even better. 

Fig. 7(a) gives (J2+3 and (Tlo for t-quark production in the range 1 < y/s/2m < 1.4 
for Jz = and = ±2, without angular cuts. Fig. 7(b) is the same except with 6^ = 
30°. We notice that the angular cuts do not make a big difference. This is because 
there is no peaked behaviour in the forward/backward directions as for bb production. As 
explained earlier, the (near) absence of resolved contributions makes the angular cuts less 
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Figure 8: Unpolarized cross sections corresponding to Fig. 7(a): (T2+3 (solid line), ulo (dashed) 
and the small P approximation (dotted). 

important experimentally as well. The most interesting feature of the HOC is that just 
above threshold, the HOC to o"(+, +) completely dominate. There is no similar behaviour 
from o"(+,— ). This suggests that the = channel is ideal for maximizing the top 
cross section not too far above threshold. At any rate, this drastic spin dependence of the 
HOC is of theoretical interest by itself and could be tested at the bb threshold as well. As 
the cross section is actually a function of only \fslm (or /3) and as((5^) (times an overall 
factor of Cq/s), the only difference would be an increase in the HOC for hb relative to 
its corresponding LO term, due to an increase in a^. In fact, the only ambiguity in the 
predictions is the choice of scale in as{Q'^). Varying in the range s/4 < < 4s, for 
^/s = 400 GeV, gives in the range 0.0878 < < 0.104 and a corresponding variation 
in the magnitude of the corrections. 

Fig. 8 gives the unpolarized cross sections corresponding to Fig. 7(a). We also plot the 
small f3 (threshold region) approximation of Ref. |^. Our results agree with this approxima- 
tion just above threshold. We see that the approximation breaks down for -yi/2m^ 1.02. 
As expected, we found that almost all of the correction comes from as, i.e. cth is almost 
negligible not too far above threshold. We found the same was true for Acts, Acth- 
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VII. CONCLUSIONS 



We have obtained complete analytical results for the production of heavy-quark pairs 
by polarized and unpolarized photons in next-to-leading order. Using these expressions, 
we computed cross sections for b- and i-quark production by photons having net spin 
Jz = 0, ±2. From the bb cross sections, we determined the background to 77 H* bb 
(standard model) coming from 77 — bb{g) (2-jet) for < A1A2 >= 0.8. The HOC to the 
J2 = channel were found to be large for s 3> 4m^. For the experimental setup considered, 
the background was safely below the Higgs signal (but still sizable) for 90 < mn ^ 150 GeV, 
even after inclusion of HOC. For t-quark production, not too far above threshold, the 
dominant contribution came from the = channel. Just above threshold, the HOC to 
this channel completely dominate. 
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APPENDIX A 



Here we present the momentum parametrizations in the frame where p4 and k are 
back-to-back. We find 



where 



Pi = 


(a;i;0,- 


■ ■ , p sinip, p cosip — U2), 




(^2;0,- 


■■,0,^2), 


k = 


{^k\ ■ ■ ■ 


, cuk sin 01 cos $2, oJk cos ^1), 


Pa = 




, — sin ^1 cos 92,—uJk cos ^1 ) , 


P3 = 




• • , |p| siu'^, |p| cos'^). 



(Al) 



Es 



s + t 
T + U 



UJ2 



S + U 
2JW2 



S2 



E. 



S2 + 2m^ 



iPl = 



2^J'S2 2-\/ 5*2 

US2 — s(t + 2m^) 



2JS'2 



cosip — 



(s + u)y 



(A2) 
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in agreement with Ref. 0. For pi, p2, ps the dots represent zeros. For k, p^ they represent 
components which depend on the remaining n — 4 angles of k. Since these components do 
not contribute to [A]|M|2^3, those angles were trivially integrated over in the phase space 

APPENDIX B 

In this appendix we list the coefficients for the various cross sections. For Adavsc/dvdw 
given in Eq. |TB|, the coefficients AAi are 

AAi = 2[-l -u/s + u^/st + m^s/tu-2m^u/t^] 

AA2 = -4[4(6M/t - 4s/m - t/T)m^/t - AsT/t^ - IQs/t + 2Au/s + s/T + Atu/sT]m^/T 

AA3 = lQ[{ls/t + ?,-?,t/u)m^/t^ + {Au^/st + 2t/s- s^/tu)m^/t + t/s + u^/st] 

AAi = A[Am^u/t^ -?,s/t + Au/s]m^/T (Bl) 

For da^se/dvdw given in Eq. |16|, the coefficients Ai are 

Ai = 2[u/t-2m'^/t + sm^/tu- Am^/tu- Am^/t'^] 

A2 = A[A{l2sT/t^ + t/T + Q)m^/uT +{s'^/uT -As^/tu-A-s^T/t^u 

- I2s/t - 2t/T - 8)mVT - 4] 

A3 = lQ[l27n^s/i^u+{s/t-lA-^t/u)m^/t^ -{"is/t + l + t/u)m^/t + u/t] 

= -A[A{2slt + tlT)m^ltu+{t/T-A)m^s'^/t^u + 2t/T-A] (B2) 

We note that 

[A]A, + {t^u) = {2mf\A]\M\lo/{Nce'e%ii'^). (B3) 
For Adahox/dvdw given in Eq. the coefficients ABi are 
ABi = AAi 

AB2 = A{s + At)m^/stu + 2{s/u + Au/s)m^/s-s'^/tu-2s/u + t/u-Au/s 
AB3 = Um'^/tu + 2m^{t - u) /tu - s'^/tu - t/u 
AB4 = [8{t- u) /s -3s/t- s/u]m^/s + 5u/s~t/s 

AB5 = AsT/tu +{t/u- l)m¥/sT^ + 4(s^ - 2t^) / st - Itu/sT + t^{s- 3t) / sTu 

ABq = 2{s/t-A)m^/tu-2m^{t-u)/tu + s^/tu + t/u (B4) 

ABj = A{u-At)m^/t'^u + 2m^{u-t)/tu + s'^/tu + t/u; ABs = {I + /Ts){t/u - 1) 

For dabox/dvdw given in Eq. ^ the coefficients Bi are 

Bl = Ai; B2 = {2m^ - s)[{2m^ -u)/st-2/u]-2m^{6m^ + t)/su 

S3 = -Am^/tu + 2m'^{u-t)/tu + s^/tu + t/u; Bi = -2m^ s/tu + 2 

B5 = 2[{2m^ + s){2s/t + t/T) -t^/T'^]/u; B^ = Sm"^ /tu + 2m'^ /u - s'^ /tu - t/u 

Bj = -Am^/tu + 2m'^{t - 3n) /tu - s^/tu - t/u\ Bg = 2{s + t^/T) /u (B5) 
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For A|M|2^3 given in Eq. and Ada-Bv/dvdw given in Eq. the coefficients Ae^ are 

Aei = -16{s/u-s/t-2)m^/u-A[s2{2 + 2s/u-t/s + u/s + 2tys^ + 8tu/s^)+2s 

- Atu/s\rr?/u + S2(4s/u - 4 - 8t^/sM - 5t/s) 

Aea = -4[2m2(2/s2M + 2/s2t+l/s2-M/s2t)+6/t-t/s2 + MVs^t]/M 

Aeg = -2[%m^{l/tu-l/s2U- s/s2t^) -2{4:s/s2 + s/t-l)m^/t-^s/u-bs/t 

- {2s^/tu - su/t + 3t + 2u + u^/t) /S2] 
Ae4 = -2{2m'^s/t + 2s + u)m^ /t; Acg = 

Agg = [32m7^i-4m2(t/M + 5 + tV5« + 5t/s + 2tVs^M + 10tVs^ + 8Ws^) 

+ Ast/u - 16t - 8t^S2/su - 5ts2/s]/2; Aej = -2m^ 

Aeg = Am^{s/tu + 2/u + l/s) - 2m^{s^/tu + 2s2/u-l+t/s) - {s^/t + s + 3s2 

+ t^/s + 3tu/s)s/u; Aeg = 8(l/n- l/s2)/t 

Aeio = 4[2(1/m + 2s/s2t - 2/s2)mVM + - 3/s2 - M/ts2] 

Aeii = 8m^(s/t + t/w + 82/3) - 2m'^[2s'^/u + 2s^/t + s + S2{2t/u + 2 + Au/t + t/s + u/s)] 

- {s"^ + t'^){t + u)/u; Aei2 = m^s; Aeia = -m^(4m^s/t + 2s + m) 

Aei4 = -m^s; Aeis = 4[2mVM - 2mVt - + ^A]/-^^ (B6) 

Aei6 = ai + ai{t^u)-8{m^ + mH/s-t)/su + 2[t^{t/u + 2)-u^{u/t + 2)]/s'^S2 

where 

ai = t[4m2(7s/t + 1 + 3t/s + 9u/s) - 2s + Qsu/t + At + u]/s2us (B7) 
For |M|2_,3 given in Eq. ^ and da-Qv/ dvdw given in Eq. the coefficients Cj are 

ei = 2[16m^(t + u)/tM+16m^(s2/u + s/t + 3) + 2m2(25s2-2t + 2s2/u) 
+ S2(2s + 5t+ 10m)]/m; 62 = -8[2m2(t + M)/s2t + 3]/tM 

63 = -2[16m''(t + u)/s2tu + 8(4 + u/t - 2s/s2)mVM + 2rr?{<ot/u + 6-u/t- s/t + 3/32) 

- 23^ ju - 33t/u - 23 - 32 + 3st/s2 " 2t^/32u]/t 

64 = 2m^{2m^/t + lf] e^ = -Am^/u 

66 = -lQm^/u + 2rr?{b3/u + 7t/u + ll) +t{A3/t- 3/u + t/u+n)] ej = 2m^ 

eg = 8m^{3/t + 2 + u/3)/3u + A(33/t + A + 232/3)m^/u + 2m^{32/u-2-t/3) 

— (s^/t + s + 3s2 + t^/s + 3tM/s)s/-u; 69 = 8(S2 + 'u)/tS2M 
= -4[8m^/M + 2m2(s/M + 2)-s7M + 4s + i(: + 3M]/ts2 

611 = lQm^[{t + uf/tu + 32/3\/3 + 8m%3-t)/u + t^/3U + 3{t + u)/3 + U^/3t\ 

+ 2m^ [u{3 + u)/t + t{t + u)/3- t32/u - 732 - 2t32/ 3 - U32/ 3] + (s^ + t^) {t + u) /u 

612 = m^(8m^/s — 4m^ — s); 613 = m^(8m^/t + 4m^ — t); 614 = 2m^; 615 = 

616 = -2[22m^{t + u)/tu + 3t/u+lA + u/t]/32 (B8) 
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APPENDIX C 



We give here the bremsstrahlung integrals, /«, appearing in Eq. (|30|). They are defined 

as 

/. = ^/ dnf,; {2mY[A]\M\l^,^Cj2[A]eJ,/s^' (CI) 

i 

(see (15) ). The fi may be exphcitly expressed as functions of 9i and 62 using the expres- 
sions in Appendix A. All the integrals here are 4-dimensional (i.e. £ = in (|27|)) and are 
determined using the general forms given in Ref. |[ 
First we list the four basic integrals, 

T + U-y ^, 2S2 , ^2 
Iq = — -'6 = — 7 — o 

S2y T+U + y S2{s + t) 

4^2 1 , xs + sis + 2S2^/SX^ 2/ , , N , 2 

h = ^^^m 5 — ^^^=, xs = 4:771 {S2S + tu) + S2S 

S2ys Jx^ Xs + s^s - 2s2Jsx^ 



Define, 



4^2 xii + st- 2^stxii - . 2( +\ , + ino\ 

hi = 7^= In ———^^=, xn = 4:771 {S2 - t) + St C2 

S2^ystXll Xii + St + 2^ystXll 



z\ = 2m s+S2S—tu, Z2 = S2U—2m s—st, z^ = m s—tu, Z4 = 2m s—tu, z^ = 2m +t 

(C3) 

We may now express the remaining integrals in terms of those listed above: 

2S2Z4 I'qZi 
-'5 — — — — + 



m2(s + t)3 {s + ty 

^9 = , o A .^3 (2^1(^2 - ^)(^ + - (zl + 2z,S2s)i2m' + ^2)} + ,/'.2 (^4 + '^^'szs) 
4b2{s + ty 4{s + t)'^ 

J _ zi Il^Zi ^ _ 8S2 f Z5 1\ 25(52 -t) 

-'10 — 7 — r~A9 ~ 7T7 — TTvi -'12 — — o — : ~ '^-'11" 



[s + ty 2(s + t)' m'^S2t \xii S2J Xiit 

^ ^ 8S2Z5 ^/ii / g2^5\ 

Xii{s2 — t)m'^t t \ Xn J 

I6S2 fl2z3 , f + sl \ 4hl r „ , .yS( S2Z^ 

hi = — -H --{S2U- SZ5-3Z5{S2-t)- [1 



stxii \xiit S2Tn? I Xiist t V Xn 

'15 = ^^{{zl + 2S2SZ,){U + T) + 2fz2{s2-t)}+^-^{{z,-u^f + 2m^sz,} 



Z2 _ 16^2(24 - 



- 7r-2 lC4j 



The integrals were put into the above form using REDUCE. The integrals not listed here 
(including the n-dimensional ones not given in Ref. ^ are straightforward and have been 
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substituted directly in (pOl). As an aside, we point out that xuit ^ u) vanishes for t> = 1/2, 
w = Wi- Hence one must avoid reaching exactly the lower bound (as for the upper) of the 
w integral, in numerical calculations. 
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Erratum: Heavy-quark production by polarized and 



unpolarized photons in next-to-leading order 
[Phys. Rev. D 51, 4808 (1995)] 

B. Kamal, Z. Merebashvili and A. P. Contogouris 



PACS number(s): 13.65.+i, 13.88.+e, 14.65.-q, 14.70.Bh 

Due to a trivial algebraic error, Eq. (30) is incorrect. On the fourth line one should 
make the replacement —1 — > 2\nx. This gives for Eq. (30) 



dvdw 



+ Sna^CpC^^l- ({2m' - s) (2 Inx (2 In ^ - - + l) 
dvdw sp \ I \ e J 



Lio 




+ 2s(3 



, sv 1 
l-21n — + - 

e 



The essential conclusions remain unchanged as the features of the figures are the same. 
Various checks were done to ensure that the analytical results (not completely presented in 
any other works) are now correct. We now agree exactly with Table 1 of Kiihn et. al. [3] and 
the corresponding polarized results in Table 1 of [G. Jikia and A. Tkabladze, Phys. Rev. D 
54, 2030 (1996)]. We also agree exactly with the unpolarized squared amplitude of [J.F. 
Gunion and Z. Kunszt, Phys. Lett. B 178, 296 (1986)] (although our 3-jet cross sections 
had a minor numerical error which made the ones presented in Figures 4 and 5 a bit too 
large on average) and all the bremsstrahlung integrals, including the n-dimensional ones, 
have been checked numerically. Finally, we agree exactly with the unpolarized virtual+soft 
cross section of [2] and consequently obtain the correct threshold behavior. The soft part of 
[2] may be obtained by making the substitution ln(st'/m^) ^ In 5 in the above expression, 
as may be easily derived. 



